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The Graphical Determination of Fresnel's Integrals. 331 

projectile if the velocities and deviations of the axes are equal. The mass of 
the shot, oc d ?; moment of inertia oc d h ; resistance x d 2 ; tilting couple 
x d 3 ; hence the velocity of precession, which varies as the tilting 
couple—moment of inertia, x d~ 2 . 

The curvature of the path due to the normal component of the resist¬ 
ance x resistance-—momentum of shot x d 2 /d 3 cc d~\ Eadius of the helix 
x curvature of the path-r-velocity of precession x d" 1 x d 2 x d. The radius 
of the helix of a 0*3-inch bullet would therefore be of the order of 
0'*4"/40 ■= 0*01 inch, which would not be detectable on the screen records.] 


On the Graphical Determination of Fresnel's Integrals. 

By John H. Shaxby, B.Sc. (Bond.), University College of South Wales and 

Monmouthshire. 

(Communicated by Principal E. H. Griffiths, F.B.S. Received March 9,— 

Read April 22, 1909.) 

The functions cos Q-?nr 2 ) dx and sin (£? tx 2 ) dx, known as Fresnel's inte- ’ 
Jo Jo 

grals, have usually been evaluated by means of converging series, for instance 
those of Cauchy and Knochenhauer. The graphical method of integration, 
by the summation of the areas contained between given ordinates, the arc of 
the curve they cut off and the #-axis, is readily applicable in this case, the 
quantity cos |-7r^ 2 (or sin |7nr 2 ) being plotted as ordinate against x as 
abscissa. The area can, then be determined between any given limits by 
Simpson's Rule. In practice, however, it is simpler and more accurate to 
apply Simpson's Rule directly to the calculated values of the ordinates, 
without plotting a curve. 

The curve y — cos \irx 2 has its zero values at x == 1, 3*, 5* etc., and is thus 
marked off into a series of loops of steadily decreasing area. The simplest 
method of computation is to find the areas of these separate loops (the first 
being only half a loop) ranging, in values of x, from 0 to 1, 1 to 3*, 
3^ to 5*, etc. The range for a particular loop will be termed its base-line. 

In the subjoined tables the calculations have been made by dividing each 
base-line into 10 parts and calculating the ordinates corresponding to the 
points of division, e.g., OT, 0*2, 0 , 3...0 , 9 for the first half-loop. 

The area of any loop may then be written as led, where d is the length of 
the base-line and h a factor depending upon the loop considered. The 
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332 Mr. J. H. Shaxby. On the Graphical [Mar. 9 ? 

values of k are shown in Table I for the first six loops and for two other 
loops of higher order. The table shows that k rapidly approaches a constant 
limiting value, 0*6366. Thus after a very few loops the area of a loop is at 
once given by the product of this constant and the length of the base-line. 
The theoretical limit of k is discussed in Section 4. 

The integration cos (-^-wx 2 ) clx to any upper limit x == n%, where n is an odd 
Jo 

integer, is thus readily performed by adding to the integral for the first few 
loops (say up to 11*) the value for the higher loops JcZd, where %d represents 
the quantity obtained by summing up the base-lines, attributing. to them 
positive and negative values alternately (because of the alternating signs of 
the loop-areas), and writing k = 0*6366. 

Table I exhibits also the areas of the first six loops, and so provides the 
data for such a calculation. 

The values of [ sin (Jw 2 ) dx can similarly be calculated for the loop-areas 
Jo 

of y = sin §ttx 2 , where the base-lines range from 0 to 2*. 2* to 2, 2 to 6*, etc., 
the values of x for which sin \irx 2 = 0. 

These loop-areas and the corresponding values of k are shown in Table XI; 
we again find the rapid approach of k to the constant value 0*6366. 


Table I.—Constants and Areas of Loops of y — cos \irx 2 . 


Limits. 

Base-line. 

Area. 

7c. 

0—1 

1 *0000 

0 *7799 

0 *7799 

1—3* 

0 *7320 

0 *4589 

0 *6268 

3*—5* 

0 -5040 

0 *3198 

0 *6344 

5*—7* 

0 *4097 | 

0 *2603 

0 *6353 

7*—3 

0 *3542 

0 *2254 

0 *6362 

3—11* 

0*3166 

j 0 *2015 

0 *6363 

99*—101* ^ 

0 *1000 

j 

0 *6365 

999*—1001* 

0 *0316 

i — 

I 0 *6366 


Table II.-—Constants and Area of Loops of y = sin |-t nr 2 . 


Limits. 

i Base-line. i 

j i 

Area. 

7c. 

0-2* 

1 

1 *4142 

0*7145 

0 *5052 

2*—2 

! 0 *5858 

0 *3706 

0 *6327 

2—6* 

0 *4495 

0 *2856 

0 *6353 

6*—8* 

j 0 *3789 

0 *2409 

0 *6357 

. 8*—10* 

0 *3339 

0 *2L24 

0 *6362 

10*—12* 

0 *3018 

0 *1921 

0*6363 
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1909.] Determination of Fresnel's Integrals . 333 

2. Integration can of course be performed to any upper limit x by applying 
Simpson’s Rule to the portion of a loop between and x, where n is the 
largest odd integer less than x 2 in the case of cos \nrx 2 , the largest even 
integer less than in the case of sin \irx 2 . A simpler method is given in 
Section 5, for cases in which n is not small. 

Table III gives the values of j cos(|-7 rx 2 )dx for x = 01, 0*2, etc., ... 1,, 

Jo 

calculated without further subdivision of the base-line of the first half loop, 
Gilbert’s values being added for comparison. It is interesting to note, with a 
view to the availability of this procedure in cases where other methods may 
not be available, that such slight subdivision can give satisfactory results. 


Table III. 


Upper limit x. 

Cx 

J COS {^TTX^dx. 

Jo 

Gilbert’s table. 

0*1 

0 *1000 

i 1 

0 -0999 ! 

0*2 

0*1999 

0 -1999 

0*3 

0 *2994 

0 *2994 

0*4 

0 -3975 

0 -3975 

0*5 

0 -4923 

0 *4923 

0*6 

0 -5811 

0 -5811 

0-7 

0 *6597 

0 -6597 

0-8 

0 *7229 

0 -7230 

0-9 

0 -7648 

0-7648 

1*0 

0 -7799 

0 -7799 


3. The values of cos (|- 7 rx 2 ) dx and sin (^irx 2 ) dx in the special case 

J x l J 

when x\ and x 2 are not small and differ but little from each other, can be 
directly found as follows:— 

d ( x 2 ) = 2 xdx and 2x — X\-\-x 2 , approximately. 


So, approximately, 



cos (J 7 rx 2 ) dx = ^ cos (Ittx 2 ) d (x 2 \ 

V “ y sh + xa V2 ’ V ; 


whence 

f* 2 2 

J cos (£ 7 rx 2 ) dx — ^ ^ + (sin —sin \irx-?). 

(1) 

Similarly, 

r. . 2 

J sin (§ 7 rx 2 ) dx — ^ ^ + ^ (cos ^irx i 2 —cos 1 7r* 2 2 ). 

(2> 
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4. Special cases. 


.. p a 

—O) COS ( J 

V x 1 


TTX 2 ) dx. 


In the case when Xi 2 = 2n — l, where n is any integer not too small for the 

2 

approximation of Section 3, the integral becomes —--^ (sin ^ttx 2 2 ±1), 

7T [Xi ~h X 2 ) 

where the positive sign refers to even values of n, the negative to odd values. 

4 


If, in addition, x 2 2 = 2n+l, the value reduces to + 


■ (%1 + X2) 


case x 2 2 —Xi 2 


2, to 


+ 


2 (X2 — X1) _ ± 0 .g 366 ^2-X!). 


,. or, as in this 


(3) 


This expression therefore gives the area of the loop of the curve comprised 
between the limits x x = (2n—l)^ and x 2 = (2tt + l)* and is numerically equal 
to 0*6366 of the base-line of the loop. 

p 2 

(i b ) sin (|? rx 2 ) dx. A similar investigation shows that if x x 2 = 2 n, the 

J x x 

9 

value becomes-~ (± 1 — cos |-7nr 2 2 )> where the positive sign is to be 

7T(^i + ^ 2 ) 

associated as before with even values of n. 

If we consider one loop, i.e. if x 2 = 2n + 2, the integral becomes again 
2 (x 2 —xi)/i r, showing that the limiting value of k is again 0*6366. 

rx r Vn p 

5 . cos (%ttx 2 ) dx = oos{%7rx 2 )dx+\ cos(|-7 rx 2 )dx, where n is the 
Jo Jo J Vn 

greatest odd integer less than x 2 . 

From equations (1) and (3) computations may now be made rapidly; thus, 
for example:— 

p*i rs p-i 

cos(-|t TX 2 )dx = I cos (- 2 7 rx 2 ) dx + I cos ( J 7 rx 2 )dx 

Jo Jo J3 

= 0-7799 - 0-4589 + 0-3198 - 0-2603 + 0-2254 + 2 ( sm f ^ . 

7r(o-0+O"l) 


The last quantity = — 0 - 0444, so the integral is 0-5615. 

rx 

sin (|- 7 rx 2 ) dx can be computed in a similar way, the integer corresponding 

Jo 

to n being even instead of odd. 

rx 

Table IV shows the values of the integral cos {^irx 2 )dx calculated by the 

Jo 

above method for upper limits ranging from 3*0 to 3*5 by steps of 0*1. 
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Table IV. 


Upper limit x. 

i 

f‘ r 

i cos (^irx^dx. 

Jo 

Gilbert’s table. 

3-0 

0 *6059 

0 *6057 

3-1 

0 *5615 

0 *5616 

3-2 

0 *4662 

0 *4663 

3*3 

0 *4058 

0 *4057 

3*4 

0 *4386 

0 *4385 

3*5 

0 *5327 

0 *5326 


6. The method is also capable of indicating the values (0*5) of the integrals 
between the limits 0 and co . Thus 


r 00 . r V(n 0 - 2 ) r® 

sin Qttx 2 ) dx — ydx + ydx, 
Jo Jo j V(«o-2) 


(4) 


where n 0 = 4£T, IS being any integer. 

== I ydx-\-k[{ — \/n{) + v / (^o +2)—+ {*— \/(%+4) + etc.} 

+ etc....] 

rV(n 0 ~2) 

= ydx—k {%“* + (wo + 4)“**+...} 


Jo' 


[S(no-2) j 

J o b 


fV(w 0 -2) r. oo 

rV(«o-2) £/oo 3 N-l 

= l ydx—s\'Zn~% -- S 
Jo oh i 


[*^26 6 

If = 28,1ST = 7, (w—2) = 26, I = 0*5622 and = 1*83. 

Jo 1 

00 

can be estimated as 2*61 by the method of approximate summation 

i 

of series by differential coefficients.* 

/* 00 

Hence ydx = 0*5622- J- x 0*6366 (2*61-1*83) = 0*500.. 

Jo 


* Boole, 4 Finite Differences,’ p. 90. 

















